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ABSTRACT 

We present deep near-infrared Keck/NIRC imaging of a recently-discovered z — 4.04 galaxy (Frye 
& Broadhurst 1998). This is lensed by the rich foreground cluster Abell 2390 [z « 0.23) into highly- 
magnified arcs 3-5" in length. Our H- and K'-h&nd NIRC imaging allows us to map the Balmer -t- 4000 A 
break amplitude. In combination with high-quality archival HST/WFPC2 data, we can spatially re- 
solve stellar populations along the arcs. The WFPC 2 images clearly reveal several bright knots, which 
correspond to sites of active star formation. However, there are considerable portions of the arcs are sig- 
nificantly redder, consistent with being observed ^ 100 Myr after star formation has ceased. Keck/LRIS 
long-slit spectroscopy along the arcs reveals that the Ly-a emission is spatially offset by w 1" from the 
rest-UV continuum regions. We show that this line emission is most probably powered by star formation 
in neighboring H II regions, and that the z = A system is unlikely to be an AGN. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: stellar content — galaxies: 
clusters: individual (Abell 2390) — (cosmology:) gravitational lensing 



1. INTRODUCTION 

Although z > 3 galaxies are being found in increas- 
ing numbers {e.g., Steidel et al. 1996a, b; Spinrad et al. 
1998; Dey et al. 1998), their study has been almost exclu- 
sively through rest-ultraviolet (UV) flux redshifted into 
the optical. As such, we merely know about the near- 
instantaneous star formation rate - the UV continuum 
is dominated by massive short-lived OB-stars, and hence 
is sensitive to star formation on the timescale of only 
~ 10 Myr. From imaging of local galaxies, morphologies in 
the UV can be very different from the appearance in the 
optical {e.g., O'Connell 1997). A comprehensive study 
of the underlying stellar populations should also utilize 
longer rest wavelengths, where the light from older and 
less massive (but more numerous) stars may dominate. 

A useful spectral feature in the study of stellar popu- 
lations is the 4000 A break, which arises from metal- line 
blanketing (predominantly Fe II) in late-type stars. Even 
if the stellar population is only moderately evolved, there 
will still be the Balmer break, dominated by an A-star pop- 
ulation. At high redshift, these rest-frame optical breaks 
are redshifted into the near-infrared. Spinrad, Dey & Gra- 
ham (1995) report a detection of an 'old' stellar continuum 
from infrared imaging of a z — 4.25 radio galaxy, and stud- 
ies of more modest redshift galaxies {e.g., LBDS 53W091 
at z — 1.55) suggest the presence of stars older than 
~ 1 Gyr, implying immense formation redshifts (Dunlop 
et al. 1997; Spinrad et al. 1997). 

At z w 4, the H (Accnt ~ 1.65 /^m) and K (Acont ~ 
2.2 /im) near-infrared (IR) pass-bands straddle the rest- 
frame 4000 A + Balmer break (Figs. 1 & 2) . For a spatially- 
resolved source, it is possible to map the underlying stellar 



populations through these colors, provided this break is 
not washed-out by a current star-burst. Here we present 
an infrared study of stellar populations in a recently- 
discovered lensed galaxy at z = 4.04 (Frye & Broad- 
hurst 1998, hereafter FB98). This system of arcs (Fig. 3) 
is among the highest-redshift gravitationally-magnified 
galaxies known (see also Trager et al. 1997; Franx et al. 
1997). These lensed sources are highly amplified but are 
not pre-selected for some AGN characteristic or extreme 
luminosity, nor are they preferentially color-selected for a 
young stellar population. We will show in §3.4 & §5.5 
that these arcs are likely to be lensed images of a nor- 
mal galaxy, affording us the opportunity to study stellar 
populations in the early Universe. Hence, we can address 
whether galaxies at the highest known redshifts are gen- 
uine 'proto-galaxies', observed in the throes of a massive 
initial star burst, or whether the star formation history is 
more extended and we are witnessing 'mature' systems. 

The structure of this paper is as follows: in § 2 we de- 
scribe the optical imaging with the Hubble Space Tele- 
scope (HST), and §3 presents spatially-resolved optical 
spectroscopy obtained with the Low Resolution Imaging 
Spectrograph (LRIS) on the 10-m Keck II telescope. Our 
near-infrared imaging of these arcs with the Near Infrared 
Camera (NIRC) on Keck I is detailed in § 4. The broad- 
band optical/near- infrared colors are used to constrain the 
dust extinction and stellar populations in § 5, in which 
we also address the lensing magnification and the intrin- 
sic luminosity of this z « 4 system. Our conclusions 
are summarized in § 6, and in Appendix A we present 
our near-IR/optical photometry of the nearby "straight 
arc" (Pello et al. 1991). Throughout, we assume a cos- 
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mology with a vanishing cosmological constant (Aq = 0), 
Ho = 50/i5okms~^ Mpc^^ and qa = 0.5 (go = 0.1) un- 
less otherwise stated. At z — 4.04, the Universe is only 
1.2/i5o^Gyr (2.1 ft-^o^ Gyr) old, a look-back time of 91% 
(87%) the age of the Universe. The proper motion dis- 
tance is cIm = 6.65 /i^o GIpc [cLm = 11-2 /i^o Gpc), and in 
the absence of magnification, 1" corresponds to a projected 
physical distance of 6.4/i5Q^kpc (ICT/igg^ kpc). 

2. LENSED ARCS AT HIGH REDSHIFT BEHIND ABELL 2390 

Deep ground-based imaging [e.g., Pello et al. 1991; 

Squires et o,l. 1996) and more recent HST studies {e.g., 
Bezccourt & Soucail 1997) of AbcU 2390 {z « 0.23) have 
revealed many giant arcs arising from gravitational Icns- 
ing of background galaxies by this massive cluster, which 
has a velocity dispersion of 1093 km s~^ (Carlberg et al. 
1996). In particular, an / = 18.2" elliptical^ close to 
the critical curve bisects an extended pair of thin, near- 
straight arcs as shown in Fig. 3. The Icnsing cluster ellipti- 
cal is at .2 = 0.22437 (galaxy #101190 in Yee et al. 1996), 
and its coordinates are given in Table 1. The arcs are 
'stretched' approximately along the major axis of the el- 
liptical, and the geometry relative to the critical curve sug- 
gests a 'cusp' arc (Pello et al. 1999). The brighter northern 
arc (/814 = 22.6™) is w 5" long, with the southern com- 
ponent (/814 = 23.1") « 3" in length. 

An arclct rcdshift survey by Bezccourt & Soucail (1997) 
detected a single emission line at 6137 A from the north- 
ern component (their arclet#20). Subsequent Keck/LRIS 
spectroscopy by FB98 determined the redshift of this 
arc to be z = 4.04 on the basis of several rest-frame 
UV interstellar absorption features {e.g., C II 1335 A, 
Si IV A 1394, 1403 A - see also §3.2), consistent with the 
emission line of Bezecourt & Soucail being Ly-a. In fact, 
Ly-a was seen in absorption rather than emission in the 
northern arc by FB98 in their discovery spectrum; the 
spectroscopic slit ran perpendicular to the axis of the 
arcs, intercepting only ^ 20% of the northern one. As 
we show in §3.3, the Ly-a emission is concentrated over 
a small portion of the arcs. The southern arc and one 
counter-arc have also been confirmed to be at this red- 
shift, with Ly-a in emission (FB98). Our own spatially- 
resolved Keck/LRIS spectroscopy (§ 3) reproduces the line 
emission from both the northern and southern arcs. The 
similarity in colors (i? — U > 3", objects 4&5 in Pello et 
al. 1991) and their positions relative to the caustic suggest 
that these arcs are lensed images of the same source (see 
§3.5). 

To imdertake a comprehensive spatially-resolved study 
of this system of arcs, we have used archival HST 
imaging in the optical (§2.1) and long-slit spectroscopy 
with Keck/LRIS (§3) to supplement our near-infrared 
Keck/NIRC images, which will be described in §4. 
Throughout, we have corrected for foreground Galactic 
extinction using the recent COBE/DIRBE & IRAS/ISSA 
dust maps of Schlegcl, Finkbcincr & Davis (1998) which 
have an optical reddening of E{B ~ V) — 0.11", equiv- 
alent to extinctions of A555 = 0.33", ^314 = 0.20", 

^Unless otherwise stated, we adopt the natural Vega-based magnitude system, where the conversion to flux density for a source of zeroth 
magnitude is given for the F555W & F814W HST/WFPC2 bands by /^{Vsbb = 0) = 3664 Jy & U{hi4 = 0) = 2466 Jy, and for the near- 
infrared filters by fu{J = 0) = 1578 Jy, fu{H = 0) = 1020 Jy & fv{K' = 0) = 688 Jy. For reference, the traditional ii"-band has a zeropoint of 
/i/ {K = 0) = 646 Jy on this system. 



We have re-reduced archival HST imaging of the 

Abell 2390 cluster from the program of Fort et al. (GO- 
5352) taken with the Wide Field/Planetary Camera #2 
(WFPC2; Trauger et al. 1994). The imaging was under- 
taken on 1994 December 10 UT in the F555W and F814W 
filters, approximating the ground-based V and I band- 
passes. The data comprised four orbits in F555W and 
five in F814W, with each orbit consisting of a single inte- 
gration of 2100s, taken at a telescope PA of -|-245°. The 
telescope was offset by an integer number of pixels between 
orbits - which unfortunately meant that 'drizzling' (Hook 
& Fruchter 1997) or other methods for recovering some of 
the WF undersampling were ineffective for this dataset. 
The data-quality files were used to create masks of unreli- 
able pixels. When the frames were registered, these pixels 
were excluded from the average. The offsets were deter- 
mined from the centroid of the star 7" from the lensing 
elliptical (Table 1). Cosmic rays were rejected using the 
crreject algorithm in IRAF . IMCOMBINE; a clipping at the 
3 a level was foimd to be satisfactory. 

For the image in each WFPC 2 filter we gen- 
erated isophotal models of the elliptical using the 
STSDAS . ISOPHOTE package, returning parameters PA = 
-^35°, f = 0.7 and re = 1'.'2 - equivalent to Re = 
5.5/i^Q^kpc at z = 0.224. We then subtracted this model 
from each image (Fig. 6). As noted by FB98, two smaller 
arclets close to the elliptical center are revealed, located 
on the radial line from the elliptical to the cluster center. 
One of these smaller arcs has been spectroscopically con- 
firmed to be at the redshift of the two main arcs (Frye et 
al. 1997). 

For the combined F555W image (8400s), the la- 
limiting magnitude reached in 1 □" was V555 = 
28.4magarcsec~^, and for the composite F814W frame 
(10500s), it was /814 = 27.6 mag arcsec"^. This is consis- 
tent with the predicted sensitivity based on the Poissonian 
counting statistics and the WF readout noise (5e^). The 
arcs appear in the WF3 chip, and are clearly resolved along 
their length. Careful comparison of transverse cross-cuts 
of these arcs with the V » 22" star 7" from the ellipti- 
cal reveals that the arcs are marginally resolved in their 
width (Fig. 4). The transverse extent of the arcs is « 0'.'3 
(FWHM), and the WF point spread function (PSF) has 
FWHMf« O'.'IS, so the width (in the absence of any tan- 
gential shear) would be 1.8 kpc (3.0 kpc) after de- 
convolution with the PSF. For this WFPC 2 imaging, the 
5(T limiting magnitudes are V555 = 26.5" and /814 = 25.7" 
using a circular aperture of diameter 0'.'5 to maximize the 
signal-to-noise ratio for the characteristic width of the arcs. 

Given that the arcs are marginally resolved in the trans- 
verse direction with HST, a naive model of a spherical 
source of diameter ~0'.'3, linearly stretched into the ob- 
served arcs 3-5" in length implies an average magnifica- 
tion of // ^10 17. where /i is the integrated flux ampli- 
fication (in this instance equal to the linear stretch); to 
first order, the radial (de)magnification with respect to 
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the cluster center will be small compared to the tangential 
stretching, although the shear stretching produced by the 
elliptical galaxy is a significant second-order effect. Al- 
though the magnification may not be constant along the 
length of the arcs (as we assume in this simple model), 
we expect the magnification to be smoothly varying as the 
images are some way from the caustic. This estimate of 
the magnification is comparable to the value of ~ 20 ob- 
tained from the lens model of FB98, who also show that 
the cluster potential is softer than isothermal'''. However, 
the existence of sub-structure (the prominent knots Nl- 
N4 and S1-S3 in the northern and southern arcs. Fig. 3) 
indicates that the source object has greater complexity. 
As both main arcs are significantly extended in the tan- 
gential direction, spatially-resolved spectroscopy along the 
arcs is possible even with ground-based seeing, and in the 
next section we describe our follow-up spectroscopy using 
Keck/LRIS with a long-slit aligned along the axis of the 
northern and southern arcs. 

3. KECK/LRIS SPATIALLY-RESOLVED SPECTROSCOPY 
3.1. LRIS Observations and Data Reduction 

We obtained the long-slit spectroscopy on the night of 
1997 July 01 UT using the Low Resolution Imaging Spec- 
trograph (LRIS; Oke et al. 1995) at the f/15 Cassegrain 
focus of the 10-m Keck II Telescope. The LRIS detec- 
tor is a Tek 2048^ CCD with 24 /xm pixels, although the 
CCD was windowed-down spatially to just the 900 pix- 
els illuminated by the 3' by 1" long-slit. The angular 
scale is 0'.'212pixel^^, and we read out the CCD in two- 
amplifier mode. The observations were obtained using the 
600 line mm~^ grating in first order blazed at 7500 A, pro- 
ducing a dispersion of 1.256 Apixel"^. The reference arc 
lamps and sky-lines have FWHM « 3.5 5 pixels (the spec- 
tral focus being best at the central wavelength), so for 
objects that fill the l"-wide slit, the velocity width of a 
spectrally unresolved line is FWHM » 220 — 325kms~^. 
The grating was tilted to sample the wavelength range 
5810 - 8630 A, corresponding to A^st = 1150 - 1710 A 
at z = 4.04. The long-slit observations were taken with 
a PA of -|-23° east of north, to cover as much of both 
arcs as possible (Fig. 6). An offset of A(RA) = -4'.'92, 
A(dec) — +1'.'90 from a nearby star (Fig. 3, with coor- 
dinates in Table 1) was performed to center the arcs in 
the slit. A total of 3600 s of on-source integration was 
obtained, and this was broken into two individual expo- 
sures each of duration 1800 s to enable more effective cos- 
mic ray rejection. The telescope was dithered by 21'.'6 
along the slit between integrations to aid in the removal 
of fringing in the red and the elimination of bad pixels. 
The observations spanned an airmass range of 1.00-1.02. 
The seeing had a FWHM = O'.'6-l'.'l over the course of 
the night. Spectrophotometric standard stars HZ 44 & 
Feige 110 (Massey et al. 1988; Massey & Gronwall 1990) 
were observed at similar airmass to determine the sensi- 
tivity function, although the night was rendered intermit- 
tently non-photometric by thin cirrus. 

Each frame first had the bias subtracted (determined 



from the overscan region appropriate to that ampli- 
fier), and was c;onverted to electrons through multiplica- 
tion by the gain (1.97e~ count"^ for the left amplifier, 
2.1 e~ count"^ for the right). A high signal-to-noise dark 
current frame (formed by averaging many dark current ex- 
posures) was then subtracted. Contemporaneous flat fields 
were obtained with a halogen lamp immediately after the 
science exposures, and these internal flats were normalized 
through division by the extracted lamp spectrum. Wave- 
length calibration was obtained from Ne+Ar reference arc 
lamps, and a fourth-order polynomial fit to the centroids of 
40 arc lines created a wavelength solution with rms resid- 
uals of 0.07 A. A small (« 1 A) wavelength zeropoint cor- 
rection was applied to the dispersion scale to compensate 
for drift in the grating angle between the arc spectra and 
the data frames; this was determined by measurement of 
the centroids of prominent sky lines. Sky subtraction was 
performed by fitting a fifth-order polynomial to each de- 
tector column (parallel to the slit), excluding from the fit 
those regions occupied by sources. 

3.2. Spectral Continuum Features and Redshift of the 
Lensed Arcs 

The two-dimensional spectrum shows continuum, spa- 
tially distinct from the light profile of the lensing el- 
liptical. This extended continuum coincides with the 
areas occupied by the arcs in the WFPC 2 imaging 
(Fig. 3). The continua from the arcs have several ab- 
sorption features, and a spectral slope fairly flat in 
longward of Arestl216A at z ~ 4.04 - Fig. 5 shows 
the one-dimensional spectral extractions. Both the 
northern and southern arcs exhibit the same absorp- 
tion features (see Table 2), which are consistent with 
O IA1302.2A/Si IIA1304.4A, Si IV AA 1393.8, 1402.8 A, 
Si II A 1526.7 A and Fe II A 1611.2 A at ^ = 4.039 ± 0.004. 
Other rest-UV absorption lines commonly observed in 
starbursts (such as C IV AA 1548.2,1550.8 A) lie close to 
prominent sky lines in the OH forest. We note that most of 
these absorption line wavelengths do not match any com- 
mon spectral feature at the redshift of the lensing galaxy 
{z = 0.22437), and so the continuum detection from the 
region of the arcs cannot be extended light from the fore- 
ground elliptical. In the continuum of the z = 4.04 arcs, 
there are also two prominent absorption features at 6450 A 
(as noted in FB98) and at 6210 A, which are attributable 
to foreground absorbers. The absorption profiles appear 
to be more consistent with the C IV AA 1548.2,1550.8 A 
absorption doublet at z = 3.225 & z = 3.010 than with 
lower-.2 Mg II AA 2795.5,2802.7 A systems. 

3.3. Spatially- Offset Line Emission Associated with the 

z = 4: Arcs 

Fig. 6 shows the area of the WFPC 2 image intercepted 
by the long-slit, displayed alongside the two-dimensional 
spectrogram (trimmed to the spectral region around Ly-a 
at z = 4.04). Our optical spectroscopy reveals two sites of 
prominent line emission, detected at high signal-to-noise 



®In a recent paper, Pello et al. (1999) use a different lens model, yielding lower magnifications equivalent to 8.5 ± 2.5 and 6 ± 2 for the 
northern and southern arcs (H3a & H3b in their notation) . This lower magnification is compatible with the observed morphology if the source 
is significantly elongated and its major axis is serendipitously near-coincident with the direction of linear stretch. In addition to these 'cusp' 
arcs, Pello et al. also report the discovery of a 'fold' arc (H5) which is located 30" away on the critical curve, but lies outside the region surveyed 
in our infrared imaging. 
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(SNR « 70). However, as is readily apparent from Fig. 6, 
only one of the four bright knots of the northern are (la- 
beled N4 in Fig. 3) has a prominent nearby emission line, 
and similarly only SI in the southern component is close 
to a region of line emission. However, neither of these 
knots detected in continuum in the WFPC 2 images is co- 
incident with the centroids of the emission lines seen in 
the long-slit spectroscopy, which lie to the north of the 
continuum in both cases (Table 1). The centroid of the 
northern line emission is 8'.'7 ± 0'.'2 from the center of the 
elliptical (resolved along the spectroscopic slit axis of PA 
23°), but the centroid of the nearest continuum knot, N4, 
is 7'.'6±(y.'l from the eUiptical in the WFPC 2 images, with 
the northern extreme of the continuum extending < (X'S 
beyond. Hence the spatial offset between the continuum 
and the center of the northern line emission is larger than 
the size of the PSFs and much greater than the uncer- 
tainty in centroiding. With its smaller magnification, this 
spatial incoherence is less pronounced in the southern arc: 
the Ly-a emission is 3.7 ± 0'.'2 from the center of the; ellip- 
tical, but the centroid of the nearest continuum knot (SI) 
is further from the elUptical (4'.'3 ± C'l). 

The spatially-offset line emission near SI has a central 
wavelength of (6134 ± 1) A, with the error dominated by 
the w 1 pixel (1.256 A) uncertainty in centering the emis- 
sion line in the slit. The wavelength is consistent with 
Ly-aA1215.7A at z = 4.0456 ± 0.001, within 400 km s"^ 
of the absorption-line redshift for the continuum. We first 
consider and reject the unlikely scenario where the line 
emission is not Ly-a arising from the z = 4.04 object re- 
sponsible for the continuum of the arcs, but instead is a 
physically distinct source. Not only would this have to be 
a chance projected alignment (to within O'.'l x (/i/10)^^ in 
the source plane), but also by coincidence the physically- 
distinct source would have to have a redshift such that its 
solitary emission line had a wavelength mimicking Ly-a at 
the redshift of the arcs to within ^ 0.001. As any line 
at a shorter rest- frame wavelength than 1216 A would be 
effectively extingiiished by the Ly-a forest absorption, an 
object at z < 4 would be required. This source would have 
to be extremely unusual in being capable of producing an 
isolated emission line with huge equivalent width: no con- 
tinuum is detected at the location of the line emission. The 
3 a lower limit on the equivalent width is M-^obs > 500 A in 
the observer's frame, derived from the noise in the F555W 
image within a ID" aperture centered at the location of 
the line emission. This line emission is extremely unlikely 
to be H/3A4861.3A at 2; = 0.262 as we do not detect 
Ha A 6562.8 A which would be > 3 times as strong (assum- 
ing case B recombination). Similarly, [O III] A 5006.8 A at 
z = 0.225 is ruled out: coincidentally, this redshift is close 
to that of the lensing elliptical, but the absence of the 
second line in the [O III] doublet makes this interpretation 
untenable - [O III] A 4958.9 A is one third as strong (Oster- 
brock 1989), well above our detection threshold. The only 
other remotely plausible interpretation for the line emis- 
sion is [O II] AA 3726. 1,3728. 9 A at z = 0.646; however, our 
spectrogram extends out to 11(3 & [O III] for this redshift, 
which are not observed. Even if the absence of these lines 
was on account of anomalous line ratios, we would still 
expect to marginally resolve the [O II] doublet (although 
any velocity smoothing greater than about 250kms~^ will 



wash out the splitting). Indeed, our spectroscopic long- 
slit did serendipitously intercept a line-emission galaxy 
at 2; = 1.129 (indicated in Fig. 6) 3" below the south- 
ern arc (Table 1), which clearly exhibits the characteristic 
doublet profile of [O II], with the redshift confirmed by 
[Ne III] A 3868. 7 A (Fig. 7). We do not see this doublet 
profile in the line emission adjacent to the arcs. Hence we 
conclude that Ly-a emission, associated with the contin- 
uum from the z = 4.04 arcs, is by far the most probable 
identification of the line. This interpretation also has the 
advantage that there exists a nearby source of UV photons 
(in the adjacent star-forming knots, N4 & SI) with which 
to generate the Ly-a recombination line; in § 5.5 we de- 
scribe a model for the transport of these photons from the 
star-forming regions to the areas observed in line emission. 



3.4. Spectral Profile of the Ly-a Emission 

The Ly-a line is clearly resolved, with a measured spec- 
tral width of 7.8 pixels FWHM for the SI knot, which was 
centered on the slit. This is equivalent to a velocity width 
of Avfwhm = 370kms~^ after deconvolution with the in- 
strumental width - really a lower limit as the source did 
not quite fill the l"-wide slit in the sub-arcsec seeing. The 
profile is asymmetric, with a pronounced red wing but a 
sharper decline in flux density on the blue side (Fig. 8). 
This appears to be a common feature in high-z starbursts 
with Ly-a in emission {e.g., Dey et al. 1998; Lowenthal et 
al. 1997), and is most likely due to an outflow of neutral 
hydrogen, where we only see the back-scattered Ly-a from 
the far side of the expanding nebula - only the photons 
on the red side of the resonant Ly-a emission-line profile 
can escape, with the blue wing being absorbed by neutral 
gas (mainly within the galaxy) . This also accounts for the 
400kms~^ difference in redshift between that determined 
from the Ly-a emission (the profile of which has been 
skewed to the red) and the redshift measurement from the 
absorption features (which may be more representative of 
the systemic redshift, although large-scale outflows in the 
ISM may result in a relative blueshift of these absorption 
lines). A spectral extraction encompassing the continuum 
and line emission regions shows that there is Ly-a absorp- 
tion blueward of the emission line (Fig. 8). The P Cygni- 
like profile is consistent with this outfiow model. Recent 
near-IR spectroscopy by Pettini et al. (1998) of the rest- 
frame optical forbidden lines and Balmer lines in 2; « 3 
galaxies also shows a > 300kms~^ redshifting of the Ly-a 
relative to the nebular emission lines. 

The velocity dispersion of cr = 160kms~^ inferred 
from the line- width of Ly-a is likely to be a lower limit, 
as the line profile is truncated. However, the veloc- 
ity width of Ly-Q: is unlikely to be representative of the 
velocity dispersion of the galaxy, as resonant scattering 
will broaden the spectral profile of the escaping photons 
(e.g., Binette et al. 1993). Spatially, the two regions 
of Ly-a emission are resolved along the slit (the ma- 
jor axis of the arcs) and extend over 1-2" FWHM. The 
flux in the line from SI is / = 37 x 10~^* ergss"^ cm~^ 
(N4 falls on the edge of the slit and is presumably sub- 
ject to significant slit losses). This is equivalent to a 
line luminosity of 0.68 x 10^^ (/x/10)~^ h'^^ ergss"^ (1.9 x 
10^2 (/x/10)-i 

^bo ^'^^^^ which in the absence of ex- 
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tinction of this line corresponds to a star formation rate® 
of 0.42 (/i/10)-i h7j Mq yr-i (1.1 (p/lO)-^ h^^f Mq yr-^). 
This recombination hne emission could be generated by 
the unobscured UV flux of about 10^ h'^ (3 x 

10" (/i/10)-i h^^f) 07-stars. 

Careful inspection of the spectrum did not reveal any 
emission lines at other wavelengths, either in the region 
of the continuum or in the areas occupied by the offset 
Ly-a emission. Table 3 gives the upper limits on the flux 
in high-ionization rest-UV lines falling within our spec- 
tral coverage, comparing with the typical line ratios from 
composite AGN spectra of QSOs (Osterbrock 1989), ra- 
dio galaxies (Stern et al. 1999) and Seyfert lis (Ferland 
& Osterbrock 1986). The non-detections of N V 1240 A, 
C IV 1549 A & He II 1640 A strongly favor the interpreta- 
tion that the Ly-o; arises from the Lyman continuum flux 
produced by OB stars, rather than the harder UV spec- 
trum of an AGN. 



3.5. Lens Geometry 

It is likely that the northern and southern arcs are im- 
ages of the same source, and the lens geometry suggests 
these have even parity (FB98). The locations where there 
is distinct Ly-a emission (at the northern-most tip of each 
arc) corroborate this. A further test comes from the Ly-a 
line centroids measured in our high-dispersion optical spec- 
troscopy; if the interpretation of these arcs as lensed im- 
ages of the same source is correct then we would expect the 
central wavelength of Ly-a to be identical; alternatively, 
two associated galaxies at the same redshift would pre- 
sumably exhibit a velocity difference of At; > lOOkms"^ 
if gravitationally bound. In fact, our slit PA was such that 
although SI was centered on the slit, N4 was off-center by 
2 pixels (Fig. 6) - the arcs are not exactly co-linear, so a 
compromise slit position was used to encompass as much 
of both arcs as possible. Allowing for this spatial offset of 
0'.'4, the line centroids of the Ly-a from N4 & SI are iden- 
tical to within our measurement error (Aw < 50kms^^) 
as illustrated in Fig. 9. The relative photometry of both 
arcs from the deepest WFPC 2 F814W image (Table 4) 
suggests that the ratio of the magnification of the north- 
ern arc to the southern is 1.6:1 (Am = 0.5™). This is 
consistent with the relative lengths of the arcs (5": 3"). 
Unfortunately, the proximity of the northern-most portion 
of the northern arc to the edge of the spectroscopic long- 
slit (Fig. 6) implies that the northern Ly-a emission-line 
region is probably subject to large slit losses, so we cannot 
use the measured line flux ratio between north and south 
as a further independent check on the relative magnifica- 
tion. 

4. A NEAR-INFRARED STUDY 

4.1. Observations with Keck/NIRC 

For our near-infrared imaging, we used the Near- 
Infrared Camera (NIRC; Matthews & Soifer 1994) at the 
f/25 forward Cassegrain port of the 10-m W.M. Keck I 

®We convert the line flux to a star formation rate using the Bruzual & Chariot (1993) spectral evolutionary models, as updated in GISSEL 
1996 (Chariot, Worthey & Bressan, 1996) with a Salpeter (1955) IMF of mass range 0.1 < Af, < 125 Mq. From this, a constant star for- 
mation rate of 1 Mq yr~^ will produce *LyC = 1-44 x 10^^ of Lyman continuum photons (Arest < 912 A). Assuming case B recombination, 
69% of these Lyman continuum photons will be reprocessed as Ly-a (Spitzer 1978, table 1). The number of Lyman continuum photons per 
unit star formation rate is a factor of ?s 2 higher (so the inferred SFR a factor of Ri 2 lower) than that obtained from the Kennicutt (1983) 
IMF. 



Telescope. NIRC is a 256^ InSb array with an angular 
scale of O'.'lSpixel^^. The observations were undertaken 
on 1997 August 25 UT in photometric conditions. Al- 
though the seeing in the near-infrared was <^ G'.'5 through- 
out the night, an intermittent drift in the telescope focus 
rendered the PSF of some frames significantly worse. 

A concern with using broad-band magnitudes to infer 
the continuum levels is possible contamination by emis- 
sion lines {e.g., Bales & Rawlings 1993). We used the K' 
filter (Acont — 2.1 /im), which at the redshift of the target 
galaxy lies red- ward of the 4000 A break, but avoids any 
possible contamination from the H/3 4861.3 A recombina- 
tion line at ^ 2.445 /im, unlike the usual ii'-band which 
has a longer-wavelength cut-off. The K'-h&nd is also a 
good approximation to the rest-frame i?-band. 

We observed Abell 2390 in both H (2280 s) and K' 
(2880 s), centering the field on the elliptical galaxy lens- 
ing the z = 4.04 arcs. The sky background was i^sky = 
13.85 mag arcsec^^ and K'^^^^ = 14.08 mag arcsec^^. The 
high background in any near-infrared observations neces- 
sitates short individual exposures to avoid saturation of 
the detector wells. We adopted exposure times of 5 s {H) 
and 10 s {K'), co-adding 12 such reads before writing out 
each data frame. Hence, for each frame we have an in- 
tegration of 60s {H) and 120s {K'). Between frames we 
offset the telescope, dithering in a 3 x 3 pattern of step size 
8". This produces an incoherent sequence of disregistered 
frames, ensuring that each pixel samples genuine sky fre- 
quently. There were 38 independent pointings in H and 
24 in K' . The NIRC field of view was oriented at position 
angle 142° east of north in order to avoid dithering the 
arcs onto brighter objects while retaining a suitable star 
in the fixed offset guider. 

Photometric calibration was achieved through observa- 
tions of the UKIRT Faint Standards FS3, FS27, FS29 
& FS32 (Casali & Hawarden 1992) over a range of air- 
masses, from which we determined the extinction coeffi- 
cients (0.06 mag airmass"^ and 0.12 mag airmass"^ for the 
H- and i^'-bands, respectively). The K' magnitudes can 
be transformed to the standard K filter by using a color 
term derived from the standard star photometry: 

K = K' - 0.136(if - K). 



4.2. Data Reduction 

The near-infrared data reduction involved the removal 
of two principal components: an additive term and a mul- 
tiplicative term. The additive term is comprised of the 
bias level and dark current of the detector, the background 
contributed by the sky, and the thermal signature of the 
telescope. The multiplicative term is the flat-field, depen- 
dent on the intrinsic gain variations of the detector pixels 
modulated by the illumination of the array. 

As NIRC does not offer a correlated double-sampling 
readout mode, we formed our own 'super-bias' frames with 
the same readout pattern as the science data, medianing 
many bias frames to reduce the readout noise. However, 



6 



Resolving the Stellar Populations in a = 4 Lensed Galaxy 



as the bias frames were taken some time away from our 
observations of the arcs, they were not appropriate for all 
pixels, particularly 'warm' pixels with variable dark cur- 
rent. Those pixels exhibiting large variability from frame- 
to-frame were subsequently identified and masked. How- 
ever, other pixels had a time-dependent dark current which 
varied on much longer timeseales, and could still be used 
with a background subtraction determined from the adja- 
cent frames taken closest in time. 

Following dark/bias subtraction, the frames were flat- 
fielded to remove the pixel-to-pixel variation in detector 
gain. We experimented with several flat-fielding tech- 
niques. Use of the median of the disregistered data stack 
to remove objects in the field, although producing frames 
which are cosmetically flat, does not provide a good repre- 
sentation of the true flat field. This is because it contains 
both the sky background (which is likely to be flat on 
the ~ 1' fleld of the detector), and the thermal emission 
of the telescope and instrument, which will not be con- 
stant across the field of view and is very signific;ant in the 
thermal-infrared /sT-band. Hence, this approach cannot 
yield a frame that is photometrically flat (i.e., adopting 
this technique to flat-field, the photometry of an object 
will be a function of the position on the array). We also 
obtained dome-fiats, formed by differencing lamp-on and 
lamp-off images of an illuminated dome spot. By taking 
the difference, the thermal signature of the telescope and 
dome spot should be removed, and all that remains is illu- 
mination of the array multiplied by the detector fiat-field. 
However, there may be non-uniform illumination of the 
dome spot, although this will vary smoothly with position 
on the array as the dome spot is very out of focus. Our 
third method of fiat-fielding involved taking a sequence 
of dithered twilight sky-flats of equal integration times as 
the sky intensity varied. We performed regression fits for 
the intensity of each pixel against the median counts in 
each frame (the sky level to a first approximation). Out- 
liers more than 3-(t from the fit were excluded, to elim- 
inate stars and cosmic rays. The slope of the regression 
line yields the normalized gain for each pixel, and the in- 
tercept is the non-flat thermal contribution (more details 
are given in Bunker 1996). This technique will yield the 
true flat-field provided the thermal signature is stable on a 
timescale of minutes - the duration of the twilight sky-fiat 
sequence. These twilight sky-flats, which were ultimately 
used in the calibration, differed from the dome-flats by a 
smooth gradient of w 4% across the field-of-view of the 
detector. The robustness of our fiat-field was confirmed 
by standard star photometry at different locations on the 
array. The dispersion in the recorded calibrated counts 
was less for the twilight sky flat than either the dome-flat 
or the median data-fiat. 

After removing the bias/dark-current and flat-flelding, 
the images were then background-subtracted. As the back- 
ground level varied on the timescale of a few tens of min- 
utes, the background for each frame was determined by 
several of the temporally adjacent images. The DIMSUM 
package'' was used to create object masks for each in- 
dividual frame. The ten frames closest in time were 
then combined, excluding the masked objects, to create a 
background image appropriate for that particular frame. 



Each frame had cosmic rays identified through unsharp- 
ened masking, which differentiated the hard edges of these 
events from the softer profile of real objects. 

We combined all the background-subtracted, flat-fielded 
frames, registering them using offsets determined from the 
centroid of a bright star close to the arcs of interest. Each 
pixel was rebinned onto a 2 x 2 grid, and the data reg- 
istered at the half-pixel level. All frames were scaled to 
zero airmass using the extinction term determined from 
the standard star photometry. Optimal inverse-variance 
weighting was employed, to maximize the signal-to-noise 
ratio. Known bad pixels and previously-identified cosmic 
rays were flagged for each individual frame and excluded 
from the co-addition. To eliminate other cosmic rays and 
artifacts which had escaped previous cuts, the two high- 
est and two lowest intensities measured for each registered 
pixel were ignored, and the inverse-variance weighted av- 
erage of the other 34 frames (in H) or 20 frames (in 
K') was taken. In the final mosaics, our 1 a sensitiv- 
ity in a in" aperture is H{1(t) = 25.46 mag arcsec"^ and 
K(la) = 24.92 mag arcsec"^ for the central region where 
all the images overlap and our noise is lowest. These 
limits are entirely consistent with Poisson statistics based 
on the sky background and the readnoise of the detector 
(30 e"). Hence, our 5cr detection limit for a point source 
is H{5a) = 23.8"* and K{5a) = 23.3"*. The PSF in our 
final combined H k. K' images had a FWHM of O'.'S. 

4.3. Analysis 

As with the WFPC 2 imaging, the surface brightness 
profile of the lensing galaxy was determined by fitting el- 
liptical isophotes of increasing radii, and this model was 
subtracted from the final NIRC image in each waveband. 

The goal of this work is to map out the spatially- resolved 
color gradients along the arcs. The coordinate transform 
from the sub-pixelated (0'.'075 pixel" ^) NIRC mosaics onto 
the intrinsically higher-resolution WFPC2 images was de- 
termined using approximately 40 compact sources com- 
mon to all frames. The resultant pixel scale is O'.'l pixel" ^. 
This process, as well as matching the different orienta- 
tions and pixel scales, also corrects for spatial distortions 
between the two detectors. 

The arcs are essentially linear. We collapse the flux 
along the short axis of the arcs in the elliptical-subtracted 
images, summing over a fixed aperture width chosen to 
maximize the signal-to- noise ratio (S/N). The extracted 
one-dimensional light profile gives the integrated flux at 
each point along the major axis of the arcs. This is some- 
what akin to the one-dimensional extraction of a two- 
dimensional spectrum, with the 'dispersion' dimension cor- 
responding to the linear distance along the arcs. 

We trained the aperture traces on the F814W im- 
age, in which the arcs had the highest S/N. We re- 
peated the photometry on corresponding aperture centers 
in the other wavebands. For the HST imaging, where the 
arcs are marginally resolved in the transverse direction 
(FWHM« 0'.'3), our extraction width was 5 pixels (O'.'S), 
« 1.5 X FWHM, which for a Gaussian profile maximizes 
the S/N. For the ground-based near- infrared images, the 
seeing in the final mosaics was « 0'.'5 and the ares are unre- 
solved in the transverse direction, so we set our extraction 



'^DIMSUM is the Deep Infrared Mosaicing Software package developed by P. Eisenhardt, M. Dickinson, A. Stanford & J. Ward. 
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width to be O'.'S. Fig. 10 shows the intensity of the arcs in 
each waveband in terms of hnear magnitudes, the flux per 
unit length along the arc. A surface brightness 'measure- 
ment' would be inappropriate for these linear, unresolved 
(or barely-resolved) systems. 

For each image, our aperture width for the photometry 
is fixed, with the trace following the centroid of the arcs. 
We convert our fixed-width aperture magnitudes to total 
linear magnitudes using aperture corrections for the frac- 
tion of the flux estimated to fall beyond our extraction 
widths. The cross-section of a bright star in our field was 
block-replicated to emulate the linear nature of the arcs. 
Analysis of the curve-of-growth of the profile perpendicu- 
lar to the major axis led us to estimate the aperture cor- 
rections to be —0.05™ in the HST images and —0.25™ in 
the NIRC images for the extraction widths adopted. 

To allow a valid comparison between the optical and 
near-infrared data with their different intrinsic spatial res- 
olutions, we boxcar-smooth the intensity maps by sum- 
ming the flux in the one-dimensional extraction over the 
width of the resolution element determined from the lower- 
resolution ground-based seeing (Of! 5). Hence, the derived 
colors (from our linear magnitudes, using fixed- width aper- 
tures) are now unaffected by differences in the resolution 
of the optical and near-infrared imaging. The maps of in- 
tensity as a function of distance along the arcs are plotted 
in Fig. 10. 

5. DISCUSSION 
5.1. Age and Reddening from Broad-Band Colors 

We wish to compare our spatially-resolved optical/near- 

infrared photometry with stellar spectral synthesis models 
of various metallicities. We also consider the effects of ex- 
tinction by dust intrinsic to the source and the absorption 
by foreground neutral gas (the line and continuum opac- 
ity of intervening absorbers). Lensing is achromatic, so 
the colors are not influenced by the fact that the system 
is magnified (once the different resolutions of the various 
images have been accounted for). 

At z = 4.04, the {H - K') color straddles the 
4000 A -|- Balmer break. The amplitude of the break is 
affected by the age, the initial mass function, the dust 
reddening, and the metallicity. It is most sensitive to age 
(Fig. 2) and relatively insensitive to dust, because of the 
small rest- wavelength baseline. We can break the age/dust 
degeneracy to some extent through use of all colors, which 
span a significant range of wavelengths. However, the 
age/metallicity degeneracy persists. 

Using ages derived from {H — K'), we can then use the 
optical/near-infrared colors as an excellent diagnostic of 
dust absorption. The F814W filter is centered at a rest 
wavelength of Arcst ~ 1600 A, and the pass-band lies long- 
ward of Ly-a. Hence it is very susceptible to intrinsic 
reddening by dust, but is immune to the effects of line 
blanketing by the intervening Ly-a forest. 

The interpretation of the (V555 — /814) rest-ultraviolet 
color is more complicated since the intrinsic stellar SED is 
extremely sensitive to dust extinction at these short rest- 
wavelengths. In addition, most of the F555W filter lies 
shortward of Ly-a from the arc;s. Hence, the F555W flux 
may have a line contribution from Ly-a (in emission or 
absorption), and the continuum will be severely depressed 



by z < 4 intervening Ly-a clouds (Madau 1995). 

5.2. Color Gradients Along the Arcs 

We detect the arcs in our near-infrared imaging over the 
entire spatial range observed in the WFPC 2 data. Fig. 10 
shows spatially-resolved colors along the tangential extent 
of both major arcs. There are statistically significant color 
gradients at the 3-5 cr level (the lower panels of Fig. 10). 
For the cases of the prominent knots seen in the WFPC 2 
imaging of the rest- ultraviolet, the colors arc compara- 
tively blue (relatively flat in /j, below Arest ~ 4000 A), 
with (H - K') K. 0.7™ in the knots. From Fig. 2, this 
is entirely consistent with an extremely young OB stellar 
population (< 20 — 30Myr). However, there are significant 
portions of the intcrknot regions that are much redder, 
with {H-K') colors of 1.4™±0.4™ between N3 and N4 and 
1.3™ ± 0.3™ between N2 and N3, with a firm lower bound 
of (H — K') > 1™. From Fig. 2, where ages are plotted 
against colors for different metallicities, the time elapsed 
since the last major episode of star formation in the ab- 
sence of dust is 40 - 200 Myr for metallicites oi Z = 0.004, 
and 40 - 150 Myr for solar metaflicity {Z = 0.020). This is 
completely insensitive to upper mass cutoffs greater than 
10 M0, as these massive stars will have left the main se- 
quence. 

A color comparison of the two highest S/N bands, the 
F814W and K' , shows that the (/814 - K') color is anti- 
correlated with the trace of the F814W flux (Fig. 10). This 
is to be expected since the F814W (A^st ~ 1600 A) traces 
ultraviolet-bright sites of recent star formation, whereas 
(/814 — K') is reddest in the older regions, and/or those 
with the most extreme dust reddening. 

In the northern arc, knots N1-N3 are not adjacent to 
the sites of Ly-a emission. However, the uncontaminated 
continuum slope of the knots {e.g., between F814W and 
i7-band) is relatively flat in /,y (Fig. 10), implying a young 
stellar population and suggestive of ongoing star forma- 
tion. Extrapolating this to shortward of Ly-a, the con- 
tinuum depression as inferred from the {V — I) color may 
be used to make a direct estimate of Da, the absorption 
due to intervening cosmological H I clouds. Outflowing 
neutral hydrogen intrinsic to the source may also play a 
significant role in absorbing the blue- wing of Ly-a (§ 3.4). 
Formally, the Da continuum break (Oke & Korycansky 
1982) is defined as 

/ ^ /.(1050- 1170A).,.N 

V /.( 1050 - 1 170 A)pred; 

{e.g., Schneider, Schmidt & Gunn 1991; Madau 1995); we 
use our best-fit SEDs (including dust reddening) to predict 
the F555W flux in the absence of H I absorption. Com- 
parison with the actual F555W flux yields a determination 
of Da = 0.5 ± 0.1, after correcting for the fraction of this 
filter which lies above the redshifted 1216 A break. 

The (V555 — /814) colors in the vicinity of knot N4, from 
which there is substantial nearby Ly-a emission, are signif- 
icantly bluer than for knots Nl— N3. This difference is con- 
sistent with the emission contributing to the F555W fiux. 
With this assumption, we estimate the equivalent width 
of the Ly-a emission at N4 to be Wrest ~ 30 A. However, 
as noted in § 3.3, the Ly-a emission extends further than 
N4, with much larger equivalent width over these regions 
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{Wrest > 100 A). The globally -averaged equivalent width 
is only Merest ~ 10 A (§5.5) typical of the Lyman break 
galaxies. 

5.3. Stellar Spectral Synthesis Models 

We adopt the spectral evolutionary models of Bruzual 
& Chariot (1993; hereafter, BC) as updated in GIS- 
SEL96 (Chariot, Worthey & Brcssan 1996). As noted by 
many authors {e.g., Spinrad et al. 1997) there are often 
disparities between different population synthesis models 
{e.g., Worthey 1994; BC; Guiderdoni & Rocca-Volmcrangc 
1987). However, we are interested in a wavelength regime 
where such differences between spectral synthesis models 
are least significant - the rest-optical. 

In fitting stellar populations there is a large number of 
free parameters (age, metallicity, reddening, star forma- 
tion history, IMF, etc.). As we have only a limited num- 
ber of independent color indices, we begin by fixing all but 
the age and the reddening, and explore three variations in 
metallicity and star formation history. 

As with the treatment of Sawicki & Ycc (1998) for the 
2; « 3 HDF star-forming galaxies (Steidel et al. 1996b; 
Lowenthal et al. 1997), we have constructed grids of BC 
models with which to fit the observed broad-band colors. 
We have conservatively adopted the Salpcter (1955) IMF 
with a large mass range (0.1 — 125 M©), compatible with 
that inferred for local starbursts (Stasinska & Leitherer 
1996.) 

We consider two mctallicitics, both solar metal abun- 
dance (for comparison with the older BC models; Z = 
0.020) and lower metal enrichment {Z = 0.004), compat- 
ible with the values derived for high-column-density ab- 
sorption systems at 2: > 2 (20% solar; c.f., Pettini et al. 
1997a), which arc probably sight-lines through galaxies. 

The final ingredient in the matrix of models is dust ex- 
tinction. Longward of Ly-a we only consider dust intrin- 
sic to the source (once Galactic extinction has been re- 
moved). There has been much debate about what is the 
most appropriate reddening law for high-redshift galaxies 
{e.g., Meurer et al. 1997; Dickinson 1998; Pettini et al. 
1997b, 1998). We adopt the extinction law of Calzetti 
(1997a, b), empirically derived from local starburst galax- 
ies, which may therefore be more appropriate to the in- 
tegrated light of dusty star-forming regions^ than those 
derived from individual stars, such as the standard LMC, 
SMC, and Galactic reddening laws. In particular, the fore- 
ground dust screen models employed with these extinction 
laws are likely to be inaccurate as the geometry of the dust 
and gas in a star forming region will be more complex. 

We have assembled a grid of SEDs, varying dust ex- 
tinction between E{B - y) = 0™ and E{B -V) ^ 0.5™ 
in increments of AE{B - V) = 0.005™, and using 144 
age steps from the BC models (from zero age to 2 Gyr, 

*The empirical reddening law for star-forming regions as determined by Calzetti (1997a, b) is parameterized in terms of the color excess 
E{B — V) for the stellar continuum of the OB stars. This is related to the color excess of the ionized gas (which would be inferred from nebular 
line emission observations) by E{B — V) = 0.44_Egas(-B — V). The differential reddening between gas and stars is due to the fact that the 
youngest stars (which provide the ionizing flux) are generally dustier than later stars (which contribute to the continuum, but no longer to the 
ionization of the gas) . The wavelength dependence of the reddening is 



the upper-limit on the age of the Universe at z = 4 for 
qo = 0.1). We then use our photometry for each location 
along the arcs to perform a maximum-likelihood fit to the 
SED at that point. In our fit, we exclude the F555W 
flux as it is prone to both line contamination from Ly-a 
and continuum-suppression below Arost = 1216 A. Using 
the magnitudes measured in the other three wavebands, 
weighted by the photometric uncertainty, we minimize the 
absolute to compute the age and reddening for the best- 
fit model. We have (n — 1) degrees of freedom for fitting 
a particular model, where n is the number of pass-bands 
used in the fit {n = 3). Using our subsequent determi- 
nation of Da (§5.2), we repeated the analysis using in 
addition a F555W flux corrected for intervening absorp- 
tion, but found this did not substantively aff'ect the best-flt 
SEDs. 

We have undertaken this analysis using different star for- 
mation histories, and place constraints based on extremes 
of the Bruzual & Chariot models, from single-event (in- 
stantaneous) starburst to a continuous star-formation rate. 
Although undoubtedly simplifled, these models bracket all 
reasonable star formation histories. Specifically, we focus 
on three scenarios: 

• a metallicity of 20% solar {Z = 0.004) and an in- 
stantaneous burst model; 

• solar metallicity {Z = 0.020) with the burst model; 

• and solar metallicity with a constant star formation 
rate. 

5.4. SEDs of Spatially Resolved Stellar Populations 

In this sub-section we discuss the fits of model SEDs 
to the broad-band photometry of several specific locations 
with extreme colors. We use this to show that non-coeval 
star formation histories are suggested for this 2 = 4 object, 
and we also place limits on the degree of dust obscuration. 

First we consider the reddest interknot regions {e.g.. 
Fig. 11, the interknot region between Nl & N2 at -3'.'96). 
Formally, a huge range of parameter space is allowed 
within the la confidence interval for all three models we 
have considered. The best fit of the single-burst mod- 
els suggests that the SEDs are dominated by a ~ 100 Myr 
population, with an extinction of E{B — y) ps 0.1™. A 
constant star formation rate model allows for an enormous 
range of ages with large extinction due to dust (as much as 
Ay ^2™). Although this model is statistically as viable 
as the single-burst ones, a constant star formation rate 
over more than 100 Myr would yield a detectable far-IR 
fiux, which is not seen by recent ISOCAM observations 
- the 3 a limit at Q.7 iim of 25y^Jy per 5" beam (Altieri 
et al. 1999; Lemonon et al. 1998) constrains the redden- 
ing to be E{B — V) < 0.7™. An even better upper limit 



k{X) = 2.656(-2.156 + 1.509/A - 0.198/A^ + 0.011/A^ 
k{X) = [(1.86 - 0.48/A)/-\ - 0.1]/A + 1.73 



where the flux attenuation for the stellar continuum is given by -Fobs(-^) 

A{V) = 4.88 E{B -V) = 2.15 Egi,s{B - V), differs from that for the Galactic law, A{V) = 3.1E{B - V). 



4.88 0.12 /im < A < 0.63 urn 

0.63 /im < A < 1.0 /im 
Fo(A) lO-O-^s^-^-^)*^^). Note that the reddening relation, 
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(Kunth et al. 1998) and studies of tlic z » 3 — 4 Lyman- 
brcak population (Stcidcl et al. 1996a, 1999), wo know that 
Ly-a emission is absent in about half the sample of 'nor- 
mal' galaxies. This is probably attributable to selective 
extinction of this resonant line after multiple scatterings 
(§5.5), with the large dispersion in the equivalent widths 
of Ly-a arising from the complex dependence on the ge- 
ometry of the ISM and the internal kinematics. 

As originally demonstrated by Bunker et al. (1998), and 
recently confirmed by the work of Pello et al. (1999), the 
broad-band colors (Table 4) integrated along the entire 
northern arc (Fig. 14) show that the overall spectrum is 
best fit by a more evolved (~ 50 Myr) population with ex- 
tinction intrinsic to the source of Ay ~ 0.3™. This implies 
that there are older stars outside the actively star-forming 
H II regions. In some areas, {H — K) f« 1.2™ (Fig. 10), 
indicating that ~ 100 Myr has elapsed since the major 
episode of star formation (Fig. 2). 

It is interesting to note that without the luxury of the 
spatially-resolved colors, the int(\s,rated flux would be most 
heavily dominated by the luminosity and the colors of the 
recent starburst regions. This would lead to the erroneous 
interpretation that the object is a primaeval galaxy, ob- 
served during its first episode of star formation. This pa- 
per shows that this is demonstrably incorrect. 

5.5. The Origin of the Line Emission 

An unusual aspect is that the Ly-a line emission is not 
coincident with the UV-continuum seen in the HST imag- 
ing (§ 3.3). The line emission regions only overlap partially 
with the northern-most edges of the two arcs, and extend 
~ 1 — 2" beyond. Therefore, the localized equivalent width 
for much of the Ly-a-bright region is huge, as the con- 
tinuum is undetected. The local 3 a lower-limit on the 
equivalent width is Wobs > 500 A in the observed frame. 
This corresponds to Wrest > 100 A (3o") in the rest-frame 
a.t z = 4.04. This lower bound is still consistent with the 
uppermost-end of the theoretical Ly-a equivalent widths 
from stellar synthesis models of star-forming regions {e.g., 
Chariot & Fall 1993) for an extremely young stellar popu- 
lation (< 10 Myr) whose UV flux is dominated by 0-stars. 
However, the observed Ly-a emission from star-forming 
galaxies is invariably much weaker when spatially aver- 
aged over the entire rest-UV extent of the galaxy. The 
Ly-a emission is typically Wrest = 5 — 30Aatz^3 {e.g., 
Steidel et al. 1996a; Warren & M0ller 1996), or even in 
absorption. Such large equivalent widths as we see from 
the emission-line regions of these lensed arcs are observed 
in AGN, although the absence of strong N V 1240 A & 
C IV 1549 A emission from our spectroscopy (§3.4) and 
the lack of X-ray emission^ ^ argue against the hard UV 

®The 3 <J flux limit in the 15" SCUBA beam (encompassing both arcs) at 850 /im of --^ 4 mjy constrains the luminosity to be 
■t-sOMm < 3 X 10" (/i/10)-i h~g Lq (<9 X 10" (/i/lO)"! Lq), following the methodology of Hughes et al. (1998). The total star forma- 
tion rate must be SFRpiR < 70 (/i/lO)"! h~g Mq yr"! (< 200 (p/lO)-'^ h^^ Mq yr"!) and the dust mass Mdust < 2 X 10^ (/i/lO)"! Mq 
(< 5 X 10^ (m/10)"^ h-^tf Mq) for /3 = 1.5, T^^^t = 50 K. 

-'"An absolute limit on dust obscuration in the high-z Universe is provided by the COBE spectrum of the microwave background, which 
is a perfect 2.73 K black-body to within 0.3% (Mather et al. 1994). Signiflcant thermal dust emission from high-z galaxies would render the 
spectrum non-Planckian {e.g., Blain & Longair 1993), and Mannucci & Beckwith (1995) show that an average extinction of only Ay > 1™ 
would be in conflict with the COBE measurements (for star formation at z > 5 and a dust temperature of T^uat = 30 K). 

^^The z = 4.04 arcs lie outside the minimum X-ray contour in the 28 ks ROSAT/HRI image of Pierre et al. (1996). The 3 a detection threshold 
within a 5" aperture (matched to the ROSAT PSF) is 6 x 10"^'' ergss"' cm~^, adopting the conversion of 7.0 x 10~" ergss"' cm~^ count"^ 
for 0.2 < E < 2 keV, Ar(H I)= 6.7 x 10^0 cm-^ and a spectral shape ax = 1 (where fu oc u "^). For a magnification of fi, the 
3 a X-ray limit constrains the source as being Lx < 0.5 (/^/10)~' on the QSO luminosity function of Boyle et al. (1993), where 



on the dust obscuration comes from the SCUBA sub-mm 
cluster lens survey of Blain et al. (1999). Comparing the 
star formation rate inferred from the rest-UV continuum 
(§ 5.6) to the upper-limit on the SFR from the SCUBA 
non-detection^ indicates that the total obscuration of the 
galaxy at I6OOA (the observed /814-band) cannot exceed a 
factor of w 10 - 15, so E{B -V) < 0.25™ for the Calzetti 
reddening law. 

Independently, Pello et al. (1999) have performed SED 
fits to the integrated flux of the arcs, and also conclude 
that for this system, continuous star formation models are 
less satisfactory than a burst scenario. For our more fa- 
vored model (passive evolution after an episode of star 
formation), the lower limit on the age is ^ 40 Myr at 95% 
confidence if we assume Ay <, 0.5™ (for both metallici- 
ties), and ^ 80 Myr in the absence of reddening (Fig. 14). 
There is significant degeneracy in the models between age 
and dust, but general conclusions can be drawn by the 
direct comparison of the results for knots and interknot 
regions within the context of each model grid. 

Comparing our grid of models with the highest surface 
brightness (knot) regions yields systematically younger 
ages than for the inter knots. Clearly, the amplitude of 
the break between the H- and K'-bands is much smaller 
in the knots (compare the SEDs of Fig. 12 with Fig. 11), 
indicating an intrinsically younger dominant stellar popu- 
lation indeed, the presence of Ly-a in the vicinity of knot 
N4 demands this. As we see Ly-a in emission nearby, we 
know the upper mass cutoff must exceed 10 Mq for there 
to be sufficient ionizing Lyman-continuum flux, which re- 
quires that star formation has occurred within the last 
10 Myr. However, in order to fit the rest-UV optical colors 
with the Calzetti reddening model, a higher extinction is 
preferred in the star-forming knots than in the older inter- 
knot regions. One could attribute this either to uncertain- 
ties in the true extinction law or, perhaps more physically, 
to a dusty environment in these star-forming regions. In 
the local Universe we know that star formation activity 
is predominantly confined to molecular clouds with high 
dust content. Recent sub-mm work on high-redshift star- 
forming galaxies (Hughes et al. 1998) also suggests a large 
extinction in the rest-UV, with re-radiation at thermal 
wavelengths by dust perhaps contributing significantly to 
the diffuse far-IR background (Schlegel et al. 1998; Hauser 
et al. 1998). Indeed, just Ay = 0.1™ of optical extinction 
win shift 25% of the rest-frame UV flux into the far-lRi°. 

Study of a prominent knot without Ly-a emission (N2) 
shows that it is still well-fit by a relatively young stellar 
population (Fig. 13), comparable to that of N4. Indeed, 
the absence of Ly-a emission is not incompatible with an 
actively star-forming region; based on local starburst work 
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spectrum of an active nucleus. The velocity width of hy-a 
is also significantly narrower than seen in the broad-line 
regions of AGN, even after correcting for the self-absorbed 
blue wing (§3.4). 

We hypothesize that the observed line emission is actu- 
ally powered by the neighboring H II regions. Lyman con- 
tinuum photons are generated in the star- forming knot, 
N4 (and its corresponding southern image, SI), consistent 
with the spectral energy distribution (SED) fits detailed 
in § 5.4. Either these Lyman continuum photons escape 
along sight-lines relatively free of gas and dust, or (per- 
haps more likely) these ionizing photons would be repro- 
cessed by H I as Ly-a, with their escape paths greatly 
lengthened through resonant scattering. Only those Ly-a 
photons which random-walk to chemically unenriched re- 
gions {i.e., areas where star formation has yet to occur) 
will escape quenching by dust. This may favor the diffu- 
sion of Ly-a photons northward of N4 & SI, away from the 
body of the galaxy. Clouds of neutral hydrogen surround- 
ing the galaxy could then act as external 'mirrors' for the 
Ly-a and any Lyman continuum photons leaking from the 
galaxy, and it is from these regions that we observe the 
line emission; if the scattering mechanism is primarily H I, 
then the cross-section of the Ly-a photons is ~ 1000 times 
greatcir than for the UV continuum, explaining why scat- 
tered continuum is not detected. Similar effects have been 
reported in radio galaxies {e.g., Villar-Martm, Binette & 
Fosbury 1996) and damped Ly-a systems (M0ller & War- 
ren 1998). 

To be more quantitative about this model, we consider 
the efficiency of creation and transport of the Lyman 
continuum photons necessary to power the observed Ly- 
a recombination line. We consider the southern Ly-a 
emission, as its measured flux is less likely to be sub- 
ject to spectroscopic slit losses than the northern region. 
The adjacent star-forming continuum knot, SI, has a 
magnitude of / 24.8™, which is a luminosity den- 
sity of LisooA ^ 7 8 X 10^7 (/z/10)-i/i^n^ erg s-i Hz- 



-i/i50^ergs 

(22 X 1028(^/10)-i/i-2^rgs-iHz-i). In the ab 
of extinction of the line. This would 



corre- 



a star 



formation 

^1 



rate^^ of SFR(Sl) 



sence 
spond to 

1.0 (a*/10)-i h^,[ Mq yr-i (2.8 (/i/lO)"! Mq yr 
and a production rate of Lyman-continuum photons 
of I-Lyc « 1-4 X 1053(^/io)-i/i^o2pi^otoj^ss-i x 

10^^ (/i/10)~^ /i^o^ photons s~^). However, assuming 
case B recombination, the observed Ly-a line flux (§ 3.4) 
requires only 40% of these Lyman continuum photons 
from SI. This could be provided by an unobscured light 
cone centered at SI, directed towards the Ly-a H I cloud, 
with an opening angle « 160° in the source plane, which 
is only « 60° in the image plane. This light cone would 
not be directed towards us, otherwise we would see Ly-a 
emission coincident with the UV continuum. Of course, 
if ionizing photons or resonantly-scattered Ly-a from the 
other (further) regions of continuum emission also reach 
the H I cloud where the Ly-a emission is observed, then 
the efficiency can be even lower - the total efficiency in the 



generation of Ly-a from the UV continuum of the entire 
southern arc (/'= 23.1'", Table 4) need only be 9%. This 
is comparable to a recent estimates by Bland-Hawthorn 
& Maloney (1999) and Dove, Shull & Ferrara (1999) of 
the fraction of Lyman-continuum photons from OB asso- 
ciations escaping the H I disk of the Milky Way. Dove 
et al. argue that the H II regions are density bounded in 
the vertical direction (perpendicular to the disk), and it 
is the escape of these photons through dynamic chimneys 
and supcrbubbles which keeps the diffuse, ionized medium 
(the Reynolds layer) ionized to scale heights of ~ Ikpc. 
It is conceivable that similar mechanisms are at work in 
the z = 4: galaxy. 

Using the spatially-integrated flux in the F814W fll- 
tcr extrapolated to the continuum just above Ly-a (the 
F555W filter is subject to Ly-a forest blanketing over 
much of its bandpass, §5.1) implies that the globally- 
averaged equivalent width of the Ly-a emission is only 
50 A in the observed frame. This value of Wrest = lOA 
is typical of the Lyman break galaxy population at high 
redshift {e.g., Steidel et al. 1996a). 

We now consider the transport of the photons. The ra- 
diation must travel ~ 1" beyond the continuum-emission 
regions of the galaxy before being scattered by H I as Ly- 
a into our line-of-sight. However, this distance is along 
the direction of magnification for these linearly-stretched 
arcs, so the physical separation is not large on account 
of this amplification 0.6 (/i/10)"^ kpc (1.1 (/i/10)"^ kpc). 
We note that without the luxury of this magnification, 
with ground-based resolution we would not be able to re- 
solve the spatially-distinct continuum and line-emission 
regions - it would look like a normal Lyman-break star 
forming galaxy, with a typical Ly-a emission line equiva- 
lent width (integrated over the entire galaxy). Even with 
HST, the sub-kpc spatial separation between the line- and 
continuum-emission would be at the very limit of the WF 
resolution (O'.'l) if the galaxy were not lensed. 

5.6. Intrinsic Properties of the Lensed Galaxy 

It has been suggested that there is a limiting surface 
brightness of star- forming regions, equivalent to a lumi- 
nosity surface density of 2 x 10^^ Lq kpc~^ (Meurer et 
al. 1997). We consider the knot N4 in the northern 
lensed arc, which is mostly likely the site of star for- 
mation responsible for the Ly-a emission. We measTuc 
the area occupied by the UV continuum of this young 
knot, rather than the area of the Ly-a emission, which 
is probably scattered significantly (§5.5). After decon- 
volution with the WFPC2 point spread function, the 
solid angle subtended by the star-forming knot N4 is 
0.05arcsec^ (out to the half-light radius), corresponding to 
0.2 (m/10)-i h^Q kpc^ (0.6 (/x/10)-i h^^ kpc^) for qo = 0.5 
{qo = 0.1), and its magnitude is K' = 22.6™. equivalent 
to Mb = -18.9™ + 2.51ogio(M/10) + dlog^oh^o {Mb = 
-20.0™ + 2.51ogio(M/10) + 51ogio/i5o) or Lb = 5.7 x 



lO9(/x/lO)-i/i5-o^L0 {Lb 



16 X 10^ {ii/lO)-^ Le), 



+5.48 (Allen 1973). The mean surface 



10«-9 (1 • 



where Mbq 

-1^ — .(.go -r z)^'* ergss"^ for go = 0, determined from QSOs out to 2 » 3. 

^^The relation between the flux density in the rest-UV around ~ 1500 A and the star formation rate (SFR in MQyr~^) is given by 
Luv = 8xl02''SFRergss-iHz-i fromMadau, Pozzetti & Dickinson (1998) for our adopted Salpeter (1955) IMF with 0.1 Mq < M* < 125 Mq 
(§ 5.3). This is comparable to the relation derived from the models of Leitherer & Heckman (1995). However, if the Scalo (1986) IMF is used, 
the inferred star formation rates are a factor of three higher for a similar mass-range. 
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brightness is therefore 1.1 x 10^^ Lq kpc~^, within the hm- 
iting value suggested by Meurer et al. (1997) but signif- 
icantly less than that derived for the star-forming knot 
in the lensed arc Gl {z = 4.92) behind CL1358-F62 
(3 X 10i2^g]jp(,-2. Franx et al. 1997). Our result is in- 
dependent of the actual magnification, as the sources are 
resolved and surface brightness is conserved, and is also 
independent of cosmology - surface brightness is subject 
only to (1 + z)"^ dimming. 

Although undoubtedly simplified, lens models (§2.1; 
FB98; Pello et al. 1999) can be used to recover the in- 
trinsic properties of the source galaxy. The inaccura- 
cies in such models are likely to be factors no more 
than a few, comparable to the uncertainties due to the 
cosmology. With the estimated linear magnification of 
fjb ~ 10, the physical scale-length of the lensed galaxy 
is r ~ 2 (/x/10)-i /i5o^ kpc (3 (/z/10)-i /i^q^ kpc), a fairly 
compact half-light radius similar to that seen in high- 
resolution HST imaging of the z w 3 Lyman break popu- 
lation (Giavalisco, Steidel & Macchetto 1996). The knots 
N4 & SI associated with the Ly-a emission occupies «10% 
of the arcs, hence we infer a size for the H II regions 
of « 200(/VlO)"i/iOT^pc (w 300(m/10)-i/i,^o^pc), com- 
parable to that observed for star-forming H II regions in 
the local Universe such as the starburst knots A & B in 
NGC 1741 (Conti, Leitherer & Vacca 1995). The dynami- 
cal crossing time for such regions is likely to be less than 
20Myr, so an episode of star formation with comparable 
duration (cohesive over these timescales) is not unphysical. 
The total luminosity density of the star-forming knots Nl- 
N4 is about £1600 A 5 ^ 1928 (/i/10)-i /i'^ ergs s'^ Hz"! 
(14 X lO^s (/i/10)"i /i5(f ergss-i Hz^i). This is equivalent 
to about 10'^ {n/l^))-'^h^^ (3x10^ (m/10)"1 h^i^) 07-stars. 
The total star formation rate inferred from the contin- 
uum around Arest = 1600 A is « 5 (/x/10)~i /igQ^ yr~i 
(14 (/i/10)^i Mq yr^i), adopting the conversion from 
Madau, Pozzetti & Dickinson (1998). This SFR esti- 
mate from the continuum is more reliable and significantly 
higher than the value derived from the total Ly-a line flux 
(§ 3.4), implying that globally this line is selectively extin- 
guished through resonant scattering by a factor of « 10 
relative to the UV continumn (§5.5). 

Introducing a dust correction to the rest-ultraviolet 
luminosity significantly increases the estimated star- 
formation rate. From § 5.4 we infer an intrinsic extinction 
of E{B -V) = 0.1™ \Av « 0.5") which corresponds to 
an attenuation in the rest-ultraviolet of a factor of three at 
1600 A. Hence, we estimate the true star formation rate of 
this stellar system to be w 15 {^/ 10)~^ h'^Q Mq yr~^ (« 
4O(/i/lO)-i/i5o^M0yr-i), which will take 3(^/10)Gyr 
(1 (ju/lO) Gyr) to assemble an L* galaxy with a mass of 
5 X 10^° /?5()^ Mq in stars, a timescale equivalent to the 
Hubble time at z = 4 and a significant fraction of the 
Hubble time at 2; = 0. Alternatively, it is likely that we 
are witnessing the formation of a sub-L* unit, which will 
in time merge with other such units to gradually assemble 
a more massive galaxy - the four bright knots may even 
be a multiple nucleus arising from a merger event. In- 
deed, Pello et al. (1999) report the discovery of the lensed 
image of another galaxy in this field at the same redshift 
(but physically distinct) 30" away. This is a projected 



distance of only 150 /igQ kpc (250 /igQ kpc) in the source 
plane, which would be crossed in « 1 Gyr at 100 km s^^. 
Hence, these two z = 4.04 galaxies may have merged by 
the current epoch. 

The absolute magnitude in the rest-ultraviolet (based 
on the /-band) is Aff|fo = -20.1™ + 2.5 logio(/x/10) -I- 
51ogio/i50 {M^em = -21-2™ + 2.5 logio(M/10) + 
51ogio/i5o)) which places it at around 0.6(/x/10) ^iigoo 
on the luminosity function determined by Dickinson (1998) 
for the z « 3 population of t/-band drop-out galax- 
ies, which have a characteristic absolute rest-UV magni- 
tude of Mi4oo* = -21™ {q = 0.5, with no dust correc- 
tion). Recent results from z » 4 Lyman break galax- 
ies show an minimal evolution in the luminosity func- 
tion from z « 3 to z « 4 (Steidel et al. 1999). From 
the rest S-band (our K' image), the absolute magni- 
tude is Mb f« -21.0™ + 2.51ogio(M/10) + 51ogio/i50 
{Mb « -22.1™ + 2.51ogio(Ai/10) + 51ogio/i5o), which 
translates to a luminosity of 3.9 x 10^° (/i/10)~i /i^g^ 
(10.8 X 10^" (Ai/10)-i/i5o^Lo). For go = 0.5, this com- 
parable to from surveys of the local Universe {e.g., 
M% K -21.0™ + 51ogio/i5o from the APM; Loveday et 
al. 1992). For q = 0.1, the i?-band luminosity of the 
lensed galaxy would be slightly above the local L*g, per- 
haps L Ki 3 (/i/10)~i i^. This almost certainly reflects 
the young age of this z w 4 galaxy, rather than a massive 
system with a mass in stars exceeding that of an L* galaxy 
at 2 = 0. 

6. CONCLUSIONS 

This paper has described a multi-waveband study of a 

z = 4 galaxy. Combining archival HST/WFPC2 imag- 
ing with deep near- infrared data taken with Keck/NIRC 
in good seeing, we have measured the spatially-resolved 
colors along a pair of z = 4.04 arcs, gravitationally lensed 
by the rich cluster Abell 2390. Comparison of the optical 
and near-IR photometry with a suite of spectral evolu- 
tionary models has enabled us to map spatially the un- 
derlying stellar populations and differential dust extinc- 
tion. Optical spectroscopy has shown that there are re- 
gions with Ly-a in emission, and that these are adjacent 
to some of the bright knots seen in the optical HST im- 
ages which sample the rest-frame UV. We have shown that 
the line emission is likely to be powered by star formation 
in these knots, rather than by an active nucleus. We see 
the Ly-Q! line morphology extending significantly beyond 
the UV continuum, which we attribute to resonant scat- 
tering from H I. In the bright knots, the SEDs are con- 
sistent with a very young stellar population (< lOMyr) 
or ongoing star formation. However, our near-IR imaging 
reveals significantly redder portions between the promi- 
nent knots, and we have shown that for these locations 
along the arcs, the most probable SEDs are those viewed 
> 100 Myr after an episode of star formation, with modest 
dust extinction of E{B — V) « 0.1™. There is degen- 
eracy in the models between dust reddening and age for 
the optical/near- infrared colors, but the most extreme sce- 
nario where the color gradients are solely due to heavy dust 
reddening of an extremely young stellar population are 
strongly ruled out by upper limits in the far-infrarcd/sub- 
mm from ISO/SCUBA. It is therefore unlikely that this 
z = 4 system in a true 'primaeval' galaxy, viewed during 
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its first major burst of star formation. Rather, our results 
suggest tliat the star formation history of this system has 
not been coeval, with current activity concentrated into 
small pockets within a larger, older structure. The evi- 
dence suggests at least two stages of star formation events 
are underway. This is not a monolithic collapse at z « 4 
- the galaxy must have begun forming at 2; > 4.5. Taken 
in the context of the recent renaissance in the study of 
normal galaxies at high redshift, our results further erode 
the paradigm of a single "formation epoch" (see also Zepf 
1997). 

Using the amplification estimated from lens models, we 
have shown that the intrinsic properties of the z = 4.04 
galaxy are comparable to the Lyman-break selected z w 
3-4 population of Steidcl et al. (1996a,1999). The current 
extinction-corrected star formation rate may be adequate 
to 'build' an L* galaxy over a Hubble time, but a more 
likely scenario may be the creation of a sub-unit which will 
undergo subsequent merging with nearby systems (such as 
the other z = 4.04 galaxy identified by Pello et al. 1999) 
to assemble hierarchically the massive galaxies of today. 
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dation. We received excellent assistance while observing 
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Goodrich, Fred Chaffee and Chuck Sorcnson. This work 
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which is operated by the Association of Universities for 
Research in Astronomy, Inc., under the NASA contract 
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with DIMSUM, and to Steve Warren for many enlight- 
ening discussions on the reduction of near-infrared data. 
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Table 1 

The coordinates of the / w 18™ lensing elliptical in Abell 2390 bisecting 

THE Z = 4.04 SYSTEM OF ARCS, AND THE LOCATIONS OF THE VARIOUS 'KNOTS' SEEN IN 
THE CONTINUUM IMAGES OF THE NORTHERN AND SOUTHERN ARCS. AlSO TABULATED 

ARE THE COORDINATES OF THE OFFSET STAR USED IN THE SPECTROSCOPY, THE 
LOCATION OF A = 1.129 [O II] A 3727 A EMISSION-LINE GALAXY WHICH WAS COVERED 
IN OUR SPECTROSCOPY, AND THE COMPONENTS OF THE NEARBY "STRAIGHT ARC." 



object 


aj20oo 




lensing elliptical 


21'' 53™ 33^52 


-M7° 4l' 57'.'4 


offset star 


21'' 53™ 33^90 


+17°4l'55'.'2 


Nl 


2l''53™33!64 


+\T 42' oo'.'e 


N2 


2l''53™33!67 


+17° 42' 01'.'5 


N3 


21'' 53™ 33!74 


+17° 42' 03'.'0 


N4 


2l''53™33!78 


+17° 42' 04'.'0 


Ly-a (N)i 
Ly-a (S)2 


2l''53™33!81 


+17° 42' 05'.'0 


2l''53™33?45 


+17° 4l' 53'.'9 


SI 


21'' 53™ 33!43 


+17° 4l' 53'.'3 


S2 


2l''53™33!41 


+17° 4l' 52'.'8 


S3 


2l''53™33?38 


+17° 4l' 52'.' 1 


II 3727 A emission galaxy^ 


21'' 53™ 33^32 


+17° 4l' 49'.'9 


straight arc A 


21'' 53™ 34!28 


+17° 4l' 54'.'7 


straight arc B 


21'' 53™ 34!41 


+17° 4l' 58'.'2 


straight arc C 


21'' 53™ 34!59 


+17° 42' 03'.'2 



^'^Exact coordinates for the line emission are uncertain as it is undetected in the imaging. 
Resolved along the spectroscopic slit axis of PA 23°, the northern line emission is 3'.'7±0'.'2 
and the southern 8'.' 7 ± 0'.'2 from the centroid of the elliptical. 

^This galaxy has two components visible in the WFPC 2 imaging; the component closest 
to the slit center is given; the other lies 0'.'4 east, 0'.'5 south. 
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Table 2 

Absorption lines in the continuum of the arcs. The mean redshift is 
z = 4.039 ± 0.04 excluding the Lv-a line. 



Aobs / A 


ID, Arost/A 


redshift / z 


6133.9 


Ly-a(cm) A 1215.67 


4.0457 


6122.1 


Ly-a(abs) A 1215.67 


4.0360 


6352.6 


Si II A 1260.4 


4.0417 


6551.9 


OIA1302.2 


4.0314 


6563.4 


Si II A 1304.4 


4.0317 


6721.0^ 


C II A 1334.5 


4.0363 


6730.6^ 


C II A 1335.7 


4.0390 


7024.7 


Si IV A 1393.8 


4.0400 


7069.2 


Si IV A 1402.8 


4.0393 


7694.0 


Si II A 1526.7 


4.0396 


7805.6^ 


C IV A 1548.2 


4.0417 


7819.2^ 


C IV A 1550.8 


4.0420 


8109. 6^ 


Fell A 1608.5 


4.0417 


8126.4 


Fell A 1611.2 


4.0437 



Affected by a nearby sky emission features. 
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Table 3 

Limits on high-ionization lines 



Emission line 


Rest-frame A 


Line fltix^ 


Ly-a flux ratio^ 


QSO ratio2 


HzRG 


Seyf 




A 


10~^^erg cm~^s~^ 


/(Ly-a)//(line) 


/(Ly-a)//(line) 


ratio^ 


ra 


Ly-a 


1215.7 


37 


1.0 


1.0 


1.0 


] 


N V 


1238.8,1242.8 


< 1.8 (3cr) 


> 20 (3cr) 


4.0 


9.7 




CIV 


1548.2,1550.8 


< 5.05(3(7) 


> 7.4^(3 cr) 


2.5 


3.9 




He II 


1640.5 


< 1.4 {3 a) 


> 26 (3ct) 


20 


5.5 





These arc 3cr upper limits in a ID" X 7A aperture. 
^From the composite broad- line region AGN spectrum in Ostcrbrock (1989). 

^From the composite high-rcdshift radio galaxy (HzRG) spectrum of Stern et al. (1999) from the MG sample. 
*From Fcrland & Ostcrbrock (1986). 

IV lies in the OH forest sky emission so the limit is less stringent. 
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Table 4 

Photometry of the arcs (integrated over 6" length for the northern arc 

AND 3" FOR THE SOUTHERN). MAGNITUDES ARE ON THE NATURAL VeGA SYSTEM, AND 
THE SIGNAL-TO-NOISE RATIO FOR EACH MEASUREMENT IS INDICATED. 



Component 


F555W 




F814W 




H 




K' 




northern arc 


24.82 (cr = 


16) 


22.61 (a = 


57) 


21.19 (cr 


= 20) 


20.38 (o- = 


25) 


southern arc 


25.14 (cr = 


17) 


23.13 (cr = 


51) 


21.99 (cr 


= 13) 


21.16 (cr = 


17) 
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Observer's frame wavelength / /xm 
12 3 



10 



8 
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4 



2 








2000 4000 6000 

Rest frame wavelength / Angstrom 



Fig. 1. — An illustration of the unreddened rest-frame optical spectra of two galaxies, one observed only 3 Myr after the end of an instantaneous 
burst of star formation (long-dash curve) and the other seen after 400 Myr have elapsed (solid line) . We also show the 3 Myr model with dust 
extinction of Av = O.S™ (using the Calzetti 1997a, b empirical reddening law), typical of high-z star-forming galaxies {e.g., Pettini et al. 1998; 
Steidel et al. 1999). Note the strong Balmer -|- 4000 A break due to the older stars. Also plotted (dotted lines) are the H and K filters in the 
rest-frame of a 2 = 4.04 galaxy, straddling the break. The SEDs come from the Bruzual & Chariot models, as updated in Chariot et al. (1996). 
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Galaxy viewed at z = 4.04 



- Z = 0.020 &c A(V) = 0.5 




100 200 300 400 500 



Age / Myr 

Fig. 2. — The evolution of the {H — K) color of a galaxy at z = 4.04 as a function of the time elapsed since an instantaneous burst of star 
formation. The solid curve is the Bruzual & Chariot model for Solar metallicity (Z = 0.020), with the dashed line showing lower metallicity, 
i solar (Z = 0.004). For this redshift, the (H — K) color is an excellent tracer of the time elapsed since the end of star formation. The dotted 
curve is the solar-metallicity model with a Calzetti-law dust reddening of Ay = 0.5"*. 
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HST / F555W HST / F814W 



o 
u 

CO 

O 




-10 





arcsec 



Keck /H (1.6 |Lim) 



Keck /K' (2.1 |Lim) 





Fig. 3.— The top panels show archival HST/WFPC 2 imaging of the cluster Abell 2390. The z = 4.04 galaxy is the arclet at PA=+23° that 
is bisected by the eUiptical. Top left is the HST V-band (F555W, 8400 s) which encompasses Ly-a, with the HST /-band (F814W, 10500 s) 
top right. The knots which arc bright in the rest-UV (and so are presumably sites of recent star formation) are indicated. Our Keck/NIRC 
images were obtained in good seeing (0.4 — 0.5" FWHM) and are shown lower left (H, 2280 s) and lower right (K', 2880 s). The radial feature 
extending from the center of the elliptical to the bottom left corner of the NIRC images is a "bleed trail" (actually an amplifier hysterisis effect) 
lying along the detector rows; we selected the PA of our observations so that this detector artifact did not ailect the photometry of the 2 = 4 
arcs. 
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arcsec 

-0.4 -0.2 0.0 0.2 0.4 




-4 -2 2 4 

pixels 



Fig. 4. — The transverse extent of the z = 4.04 arcs (diamond symbols), averaging 0'.'6 along the northern arc around N3. A Gaussian 
fit to the distribution is shown, with u = 0'.'123 (FWHM=0"3). Overplotted is an appropriately-scaled cross-cut of the nearby bright star 
(it ~ 0''078, FWHM=0'.'18) which is undersampled by the WFPC2 pixels (0''1). Therefore, the arcs are marginally resolved in their width. 
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Fig. 5. — One-dimensional spectral extraction of the continuum in the southern arc, with the wavelengths of common rest-UV absorption 
features marked. The extraction width was 8 pixels spatially (1'.'7) and the spectrum has been smooth with a box-car of 7 pixels (9 A). The 
dotted line shows the 1 a noise level, and the lower panel displays the sky spectrum. The prominent emission line in the background spectrum 
at 6711 A does not correspond to a known sky line, and is probably attributable to an LED light somewhere in the detector's light path. 
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iiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii 




Arcsec Wavelength Arcsec 



Fig. 6. Left: the F814W image with eUiptical galaxy model subtraeted (note the counter arcs perpendicular to the axis of the main 
arcs, predicted by the lens model of FB98). The area covered by the long-slit optical spectroscopy is shown (slit axis is vertical). The right 
panel is this elliptical-subtracted image, smoothed to 0''6 seeing. The apparent slight over-subtraction of the elliptical at the very center 
arises from the finite-pixel sampling skewing the central isophotes; it does not affect the photometry at the location of the arcs. Center: the 
LRIS spectrum, with the long-slit aligned along the arcs. The dispersion axis is horizontal, with wavelength increasing from left to right. The 
image has been smoothed by convolving with a Gaussian kernel of cr = 1 pixel. The positions of N4 (— 7'.'6) and SI (4'.'3) are indicated by the 
horizontal bars. Note the spatial range of Ly-o emission, which extends well beyond the detectable continuum of the arcs. The knot N4 (top) 
lies on the edge of the slit, hence the slight spatial offset of the Ly-a line centroid and the lower flux compared to the line from SI (bottom). 
Also indicated on the left panel is the z = 1.129 [O II] A 3727 A galaxy also falling on our spectroscopic long-slit (see Fig. 7). 



24 



Resolving the Stellar Populations in a ^; = 4 Lensed Galaxy 



^(rest) / ^ 

3720 3725 3730 3735 



8 



° 2 







-2 



"1 I I I I I I I I I I I I I I I I I I I r 



[OII]3727 galaxy @z=1.129 



J I I L 



J I I L 



J I I L 



J I I L 



7920 7930 7940 7950 

^(obs) / ^ 



Fig. 7. — One-dimensional spectral extraction of the [O II] AA 3726.1,3728.9 A galaxy at ^ = 1.1293, 3" below the southern arc and intercepted 
by our spectroscopic long-slit (Fig. 6). The extraction width is 7pixels (I'.'S). Note that the emission-line doublet is clearly resolved. We do 
not see this structure in the emission lines from the axes, implying that their origin is not [O II] 3727 A at z = 0.64. 
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^(obs) / ^ 



6160 



Fig. 8. — One-dimensional spectral extraction of the southern arc, showing the region around Ly-a at z = 4.04. The extraction width is 
8 pixels (I'.'V), and encompasses both line- and continuum-emission regions. The asymmetric emission line profile is readily apparent, with the 
sharp decline on the blue side due to absorption by H I within the galaxy - a blueshifted Ly-a absorption trough is visible from the outflowing 
HI). 
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Fig. 9. — The line profiles of Ly-a emission from the star-forming knots N4 & SI in the northern and southern arcs. The ccntroids are 
identical to within our measurement error (A^' < SOkms^"'), supporting the hypothesis that these are lensed images of the same source. For 
comparison purposes, the Ly-a flux from N4 (which was close to the edge of the spectroscopic long slit) is normalized to the peak of the SI 
emission, which is well-centered. The integrated line flux from SI is / = 3.7 x 10"^'^ ergscm"^ s~^. 
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Arcsec relative to Elliptical 



Fig. 10. — Magnitudes (top panel) and colors (middle panel) of the arcs. The northern arc is displayed left, with the southern arc on the right 
with the same relative scale. The (/814 — K') color has been offset by +1™ for clarity. Linear magnitudes are shown (the brightness per unit 
length), and the errors bars correspond to the la Poissonian noise when averaging over a resolution element (si 0'.'5). Note the similarity in 
the flux-per-unit length in the northern and southern arcs - these are probably images of the same source, and surface brightness is conserved 
in lensing. The bottom panel shows the statistical significance of these color gradients: for each resolution element, the number of standard 
deviations away from the mean intergrated color of the arc is displayed. The symbols are the same as the middle panel - {H — K') open 
diamonds, (/ — K') filled circles and (V555 — /814) open triangles. The variation in color is significant at the 3-5 o" level. The large deviation 
in (^555 ~ ^814) from the average at < — 7'.'7 is because of line contamination from the site of Ly-a emission. 
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Fig. 11. — The parameter space of age and dust fitting the spectral energy distribution (SED) of a portion of the hnear z = 4.04 arcs, in 
this case the interknot region between Nl & N2. The la, 2 a and 3 a contours are shown, based on the absolute distribution for n — 1 = 2 
degrees of freedom. These correspond to 68%, 95% and 99% confidence intervals. The formal best fit (lowest x^) is denoted by a circle, and the 
inset panel shows the photometry in the four wavebands plotted over the appropriately-normalized SED for these parameters. Only the F814W 
('/-band'), H & K' are used in the SED fitting; the F555W ('V-band') lies shortward of Ly-a and is subject to absorption by intervening H I. 
The fits are shown for a metallicity of 20% solar (Z = 0.004) and an instantaneous burst model (a); solar metallicity {Z = 0.020) with the 
burst model (b); and solar metallicity with a constant star formation rate (c). 
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Fig. 12. — The age and dust fits to the SED of knot N4 in the northern arc, which is a site of Ly-a emission. Refer to the caption of Fig. 11 
for details. 
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Fig. 13. — The age and dust fits to the SED of knot N2 in the northern arc, which has no detectable Ly-a emission. Refer to the caption of 
Fig. 11 for details. 
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Fig. 14. — The broad-band optical/near-IR flux from the entire northern arc. Also plotted are reddened instantaneous-burst stellar population 
models viewed at various ages, arbitrarily normalized to the flux measured from the HST/WFPC 2 F814W image. The flux in F555W is severely 
attenuated by the opacity of the intervening Ly-o forest. The colors are best reproduced by a stellar population ~ 50 Myr old, with in situ 
dust reddening of E{B — V) ~ O.!"*. Note that at 2; = 4.04, the strong Balmer + 4000A break due to the older stars lies between the H- and 
K'-Rlters. The AB magnitude system (right axis) is defined as: AB = —2.5 logjg f,, — 48.57 (Oke & Gunn 1983) where the flux density, f^, is 
in units of ergscm~-^s~^ Hz~^. 
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APPENDIX 

APPENDIX A: THE "STRAIGHT ARC" 

Our near-infrared imaging covers tlie nearby "straight arc" (Pello et al. 1991), which has generated much interest given 
the challenge to lens models in reproducing this geometry (e.g., Kassiola, Kovner & Blandford 1992). However, our H and 
K' images indicate that the "straight arc" is in fact a serendipitous alignment of two or three significantly less straight 
components (Fig. Al). The photometry of these components is detailed in Table 5. Although the colors of A & B are 
similar, C differs significantly it is redder by a magnitude in (Isu — K') - adding weight to the interpretation that these 
are physically distinct objects. Indeed, spectroscopy by FB98 indicates that the two components have different redshifts 
{ZA = 1.033 k zc = 0.913). 
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Fig. A1. — The HST/optical and Keck/near-infrared imaging of the "straight arc" behind Abell 2390. The H and K' imaging clearly shows 
that this is composed of two or three physically distinct objects, one of which (C) exhibits much redder colors than A or B. 
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Table A1 

The colors on the Vega system of the components of the "straight arc" 

LENSED BY AbELL 2390. THE SIGNAL-TO-NOISE RATIO IS TYPICALLY > 100 FOR EACH 

component in EACH BAND, SO THE RANDOM ERRORS ARE ^ 0.01™, AND THE 
SYSTEMATIC UNCERTAINTY IN THE PHOTOMETRIC ZEROPOINTS DOMINATES. 



Component K' {V555 - K') {hu-K') {H - K') 



A 19.23 1.67 0.67 0.59 

B 19.06 1.79 0.78 0.55 

C 18.00 2.77 1.72 0.92 
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